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1. Introduction 

Computation of high-speed hydrodynamics requires high-fidelity simulation to resolve flow features involving 
transient flow, cavitation, tip vortex and multiple scales of unsteady fluctuations. One example of this type in 
aerospace is related to liquid-fueled rocket turbopump. Rocket turbopumps operate under severe conditions at very 
high rotational speeds typically at thousands of rpm. For example, the Shuttle orbiter low-pressure-fuel-turbopump 
creates transient flow features associated with reverse flows, tip clearance effects, secondary flows, vortex 
shedding, junction flows, and cavitation eff ects. Flow unsteadiness originating from the orbiter Low-Pressure-Fuel- 
Turbopump (LPFTP) inducer is one of the major contributors to the high frequency cyclic loading that results in 
high cycle fatigue damage to the flow liners just upstream of the LPFTP. The reverse flow generated at the tip of 
the inducer blades travels upstream and interacts withthe bellows cavity. Simulation procedure forThlsType High- 
speed hydrodynamic problems requires a method for quantifying multi-scale and multi-phase flow as well as an 
efficient high-end computing strategy. The current paper presents a high-fidelity computational procedure for 
unsteady hydrodynamic problems using a high-speed liquid-fueled rocket turbopump. 


2. Computational Procedure 

The incompressible Navier-Stokes flow solver based on the artificial compressibility method was used to compute 
the flow of liquid Itydrogen in several test geometries including feed line and an inducer for pumping. All 
computations included tip leakage effects with a radial tip clearance of 0.006 inches, a pump operating condition of 
amass flow rate of 154.7 Ibm/sec, and a rotational speed of 15,761 rpm. One of the computational models is 
illustrated in Figure 1 which shows the LPFTP inducer with 4 long and 4 short blades, and a straight duct, which 
extends 4 duct diameters upstream of the inducer. The bull nose of the inducer, and split blades are included in the 
model. The objective of studying this model is to compare unsteady pressure values against experimental data. To 
resolve the complex geometry in relative motion, an overset grid approach is employed. The geometrically complex 
body is decomposed into a number of simple grid components. Connectivity between neighboring grids is 
established by interpolation at each grid outer boundaries. Addition of new components to the system and simulation 
of arbitrary relative motion between multiple bodies are achieved by establishing new connectivity without 
disturbing the existing grids. This computational grid has 57 overset zones with 26.1 Million grid points. 



Figure 1, Surface grids for LPFTP and the LH2 feed 
line including circumferential holes upstream of the 
pump. 



Figure 2. INS3D-MLP and MPI+OMP performance 
on Columbia platform. 



The present computations are performed utilizing the INS3D computer code 1 " 3 , which solves the incompressible 
Navier-Stokes equations for both steady-state and unsteady flows. To obtain time-accurate solutions, the equations 
are iterated to convergence in pseudo-time for each physical time step until the divergence of the velocity field has 
been reduced below a specified tolerance value. The total number of sub-iterations requited varies depending on the 
problem, time step size and the artificial compressibility parameter used. Typically the number ranges from 10 to 30 
sub-iterations. 

The computational requirement for the model problem is big, and thus two distinct parallel processing paradigms 
have been implemented. These include the Multi-Level Parallel (MLP) and the MPI/OpenMP hybrid parallel 
programming models. Both models contain coarse and fine grain parallelism. Coarse grain parallelism is achieved 
through a UNIX fork in MLP and through explicit message passing in tire MPI/OpenMP hybrid code. Fine grain 
parallelism is achieved using OpenMP compiler directives in both the MLP and MPI/OpenMP hybrid codes. 
Computations were performed to compare the scalability between the MLP and MPI/OpenMP hybrid versions of the 
INS3D code on the Columbia system using the BX2b processors. Figure 2 displays the time per iteration (in 
minutes) versus the number of CPUs for both the MLP and MPI/OpenMP hybrid versions of the code. 


3. Computed Example 


The Reynolds number for the calculations is 36 Million. The computed results show a significant time-periodic 
backrjfio w^generated by the inducer reaching 15-20% of the tip velocity and a jet flow of 10-15% of die inducer tip 
speed which penetrates into the bellows cavity creating an unsteady recirculation region in the cavity. The reverse 
flow and unsteady recirculation regions create an unsteady interaction between the duct of the feed line and the 
cavity region just upstream of the inducer resulting in high frequency cycle loading on the feed line structure. The 
back-flow also creates swirl in the cavity on the order of 10% of die inducer tip velocity. The LPFTP inducer with 
flowliner geometry is analyzed to represent the time-periodic nature of the flow field. The computed results show, 
dominant 4N unsteadiness at a fixed location as shown in figure 3(a). Figure 3(b) shows maximum and minimum 
non-dimensional pressure values recorded from die experimental data. Comparisons between CFD results and hot 
fire test data also show good correlation in the non-dimensional pressure amplitudes. 
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(a) (b) 

Figure 3. Time history of non-dimensional pressure during one inducer rotation (model 1, 14 th inducer rotation), and 
Min/Max values of non-dimensional pressure in hot fire test. 

In the full-length version, further details of the high-fidelity computational procedures will be discussed using the 
example described above. 

References 

1. Kiris, C., and Kwak, D., ‘Parallel Unsteady Turbopump Simulations for Reusable Launch vehicle,’ Frontiers of 
Computational Fluid Dynamics 2002, Caughey, D.A. and Hafez, M., ed. World Scientific, 2002. 

2. Kiris, C., Kwak, D., and Rogers, S., ‘ Incompressible Navier-Stokes Solvers in Primitive Variables and Their 
Applications to Steady and Unsteady Flow Simulations,’ Numerical Simulations of Incompressible Flows, Hafez, 
M. , ed. World Scientific, 2003 . 

3. Rogers, S. E., Kwak, D. and Kiris, C., “Numerical Solution of the Incompressible Navier-Stokes Equations for 
Steady and Time-Dependent Problems,” AIAA Journal, Vol. 29, No, 4, pp. 603-610, 1991. 





























Flowiiner Analysis 













ay 


oi 


a 


”0 nT 




o .m 


D-. O 




£ < 


■■s is 


CO erf 


O 42 


w ^ 


£3 


13 


{utui} uo!}GJ9}j jsd auuix 






